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Abstract 

We consider random Hermitian matrices made of complex or real M x N 
rectangular blocks, where the blocks are drawn from various ensembles. These 
matrices have N pairs of opposite real nonvanishing eigenvalues, as well as 
M — N zero eigenvalues (for M > N .) These zero eigenvalues are "kinemat- 
ical" in the sense that they are independent of randomness. We study the 
eigenvalue distribution of these matrices to leading order in the large N, M 
limit, in which the "rectangularity" r = 4£ is held fixed. We apply a variety of 
methods in our study. We study Gaussian ensembles by a simple diagrammatic 
method, by the Dyson gas approach, and by a generalization of the Kazakov 
method. These methods make use of the invariance of such ensembles under 
the action of symmetry groups. The more complicated Wigner ensemble, which 
does not enjoy such symmetry properties, is studied by large ./V renormalization 
techniques. In addition to the kinematical 5- function spike in the eigenvalue 
density which corresponds to zero eigenvalues, we find for both types of en- 
sembles that if \r — 1| is held fixed as N — > oo, the N non-zero eigenvalues 
give rise to two separated lobes that are located symmetrically with respect to 
the origin. This separation arises because the non-zero eigenvalues are repelled 
macroscopically from the origin. Finally, we study the oscillatory behavior of 
the eigenvalue distribution near the endpoints of the lobes, a behavior governed 
by Airy functions. As r — > 1 the lobes come closer, and the Airy oscillatory be- 
havior near the endpoints that are close to zero breaks down. We interpret this 
breakdown as a signal that r — * 1 drives a cross over to the oscillation governed 
by Bessel functions near the origin for matrices made of square blocks. 



PACS numbers: ll.10.Lm, 11.15.Pg, ll.10.Kk, 71.27.+a 



1 Introduction 



In random matrix theory, a number of authors |], 0, || £|, || || [7|, § have studied the 
eigenvalue distribution of a Hermitian matrix H of the form 

MS?). 

in which C is an iV x iV complex random matrix taken from an ensemble with the 
probability distribution 

P(C) = ^exp(-iVTrC t C) , (1.2) 

Zj 

with iV tending to infinity. These so-called chiral matrices appear in a variety of 
physical problems. For example, in quantum chromodynamics one typically inte- 
grates over the quarks and studies the so-called fermion determinant. The gluon 
fluctuations are then often treated approximately by saying that they effectively ren- 
der the relevant matrix in the determinant random [1], |9], |I0|| . The chiral structure 
corresponds to left and right handed quarks. As another example, Hikami, Shirai, 
and Wegner |TT], |T2, |TJ[ have proposed a model for electron scattering off impurities 



in quantum Hall fluids in the spin-degenerate limit. The blocks in (1.1) correspond to 
spin up and spin down electrons. In the same spirit, one may consider any problem 
involving random scattering between two groups of states, for example, between two 
cavities. As pointed out by Nagao and Slevin [|J, these matrices also appear in the 
study of transport in disordered conductors. In this paper, we study a slight gener- 
alization of this problem, with C taken to be an M x N rectangular matrix, with M 
and N both tending to infinity. For M — N of order N°, we expect the density of 
eigenvalues to be the same as for the M = N case. Here we would like to study the 
case where the measure of rectangularity, 

r = M/N , (1.3) 

is held fixed as both M and N tend to infinity. Some aspects of this problem have 
been studied before and we will note the appropriate references below. We denote 
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the matrix elements of C by 

C ia , where i — 1,2, M and a = 1,2, N. 

With no loss of generality we assume throughout this paper that M > N, namely, 
that r > 1. Our notation is such that Latin indices always run from 1 through M, 
whereas indices denoted by Greek letters run from 1 to iV.[] As a result of their specific 
structure these matrices have N pairs of opposite real nonvanishing eigenvalues, as 
well as M — N zero eigenvalues. These zero eigenvalues are "kinematical" in the sense 
that they are independent of the probability distribution. 

We derive the eigenvalue distribution of these matrices to leading order in the large 
N, M approximation for various ensembles of random blocks. We consider random 
Hermitian matrices made of complex or real M x N rectangular blocks, where the 
blocks are drawn either from ensembles symmetric under some group action or from 
non-symmetric ensembles. For concreteness, we specialize to Gaussian ensembles in 
the first case. In the second case we analyze matrices of the "Wigner Class" , namely, 
blocks whose entries are drawn independently one of the other from the probability 
distribution. We find, not surprisingly, that to leading order in the large N, M ap- 
proximation, all the ensembles we studied result in the same eigenvalue distribution. 
In addition to the kinematical ^-function spike in the eigenvalue density which corre- 
sponds to zero eigenvalues, we find that if |r — 1| does not scale to zero as iV — > oo, 
the N non-zero eigenvalues give rise to two well separated lobes that are located 
symmetrically with respect to the origin. For random Hermitian matrices that are 
not made of blocks, the qualitative universality of the Wigner semicircular eigenvalue 
distribution is well understood as a result of the competition between level repulsion 
and the fact that very large eigenvalues are suppressed. Similar arguments explain 
the universality of the eigenvalue distribution we observe here for matrices made of 
rectangular blocks. Each lobe arises qualitatively for the same reasons that lead to the 
semicircular distribution. In addition, separation of the two lobes arises because the 

1 We shall deviate slightly from this convention only in section 2 where [i, v will run over all 
possible M + N values. No confusion should arise. 
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non-zero eigenvalues are repelled from the origin by the macroscopic number (M — N) 
of zero eigenvalues. 

In this paper C'C and CC* are Hermitian non-negative matrices of dimensions 
N x N and M x M, respectively. We are interested in the expectation value of the 
resolvent 

A straightforward calculation then yields a simple relation between Gn^m(z) and the 
resolvents of C^C and CCt , 



Gn,m(z) 



N + M 



Tr,, n -= — — — - + Tr, 



z 



;i.5) 



where the subscript on each trace indicates the dimension of the matrix which is 
being traced over. The z 2 dependence of the resolvents in (|1.5|) arises because the 



eigenvalues of H in fll.l| ) occur in real opposite pairs. Indeed, given an N dimensional 
vector x and an M dimensional vector y such that ^ X ^ is an eigenvector of H for an 

eigenvalue A, then ^ X ^ is an eigenvector for —A. In other words the matrix H (the 
"Dirac" operator, with its "chiral" components C and C') anti-commutes with the 
"75" matrix ( ^ ^ Y The cyclic property of the trace implies the basic relation 
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Tr = Tr 1 . (1.6) 

(M) z 2 - CCt (iV) z 2 - CtC z 2 K ] 

This relation reflects the fact that C'C and CC^ share the same strictly positive 
eigenvalues, but the M x M matrix CCt has additional M — N zero eigenvalues. 



Combining (|1.5|) and ( |1.6|) we therefore arrive at the two alternative expressions 

r ( \ ( M — N \ 1 2z 1 
G NM {z) = - ■ — - + — ■ „ r Tr 



N + M J z M + N (JV) z 2 -CW 
N - M\ 1 2z 1 

WTm) z + MTiv r(M) ^ 2 - cct ' ( • ^ 

that allow us to express Gnm(z) solely in terms of either C'C or in terms of CC'. 
For later use we introduce the following notation 



in terms of which we rewrite ( |1 . 7| ) as 



A , , fM-N\ 1 / 2N , 



\N + MJ z \M + N 
fN-M\ 1 / 2M , - 



V.V • MJ z \M + N 

Throughout this paper G stands for an unaveraged resolvent. The corresponding 
averaged quantity will be denoted simply by G. 

This paper is organized as follows. We will first apply a variety of methods to study 
the density of eigenvalues. In Section 2 we derive the density of states of matrices H 
whose rectangular blocks are drawn either from the unitary or from the orthogonal 
Gaussian ensemble, employing diagrammatic techniques. Section 3 is devoted to 
blocks with independent entries (which we refered to |TJ[ as the "Wigner Class".) 
This ensemble is more difficult to handle, because of lack of symmetry. We overcome 
this difficulty by applying recursive manipulations of the large N renormalization 
group |15], [H| |T7|]. We find that as far as the density of states is concerned, this 
ensemble falls (in the planar limit) into the same universality class as the symmetric 
ensembles. In Appendix A we provide a proof of the central limit theorem by means 
of the large N renormalization group, as yet another example of its usefulness. In 
Section 4 we present the Dyson gas approach to these issues. After completing our 
work we realized that the results we obtained following the Dyson gas approach 



already appeared in [18|. Nevertheless, we include this section here for the paper to 
be self-contained and also because Section 5 partly relies on it. In Section 5 we first 
generalize Kazakov's method [|l9j to rederive the results of Section 2, and then use 
this method to determine the oscillatory fine structure of the eigenvalue density in 
Section 2, close to its support endpoints. We find that this oscillatory behavior is 
governed by Airy functions. As r — > 1 the lobes come closer, and the Airy oscillatory 
behavior near the endpoints that are close to zero breaks down. We interpret this 
breakdown as a signal that in the limit r — > 1 drives a cross over to the oscillation 
near the origin in the density of eigenvalues of matrices made of square blocks, an 
oscillation governed by Bessel functions. 
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2 A diagrammatic approach 



As a simple warm up exercise, and in order to set the stage, let us first apply the 
by-now well-known diagrammatic method to derive the Green's function 

in the large N, M limit. To this end, let us consider the averaged resolvent 



where the indices \i and v run over all possible M + N values. The average in ( 2.2 
is performed with respect to the Gaussian measure 

P{C) = — exp l-y/NM m 2 Tr C ] C] , (2.3) 

Zj 

where 

. M N 

Z = Yl J[ d Re C ia d Im C ia exp [-VWm m 2 Tr C ] C] (2.4) 

J i=Xa=l 

is the partition function. We have introduced a normalization factor of yMN in (|2.3|) 
so as to be consistent with ( |1.2| ) in the N = M case. This factor renders ( p.3| ) and ( |2.4| ) 
manifestly symmetric under M ^ N. Some other normalizations, not symmetrical 
under M <-> A/", can always be introduced by multiplying the parameter m 2 by an 
appropriate factor of r — M?. Borrowing some terminology of gauge field theory we 
may consider C, C< as "gluons" (in zero space-time dimensions), and G^(z) as the 
propagator of "quarks" (with complex mass z) which couple to these "gluons". We 
now proceed to calculate G^(z) diagrammatically. The two-point correlator associated 
with ( |2.3| ) is clearly 

{C ia C* p ) = (2.5) 

This expression is the gluon propagator. The bare quark propagator is simply -. The 
quark-quark-gluon vertex factor is 1. These Feynman rules are summarized in Fig. 

(!)• 
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(a) 



a p 
(b) 



i a 



a 



(c) 



Fig. (1) - The Feynman rules: (a) The bare quark propagator, (b) The gluon propagtor, (c) The 
bare quark-quark-gluon vertex. 



The weight (|2.3| ) is Gaussian, so there are no gluonic self interactions. Due to the 
block structure of H, the quark-quark- gluon vertex converts the color type carried 
by a quark. Namely, it converts a quark qi into a quark q a (and vice- versa), but not 
into a quark qj. Recall that the dominating diagrams in the large N, M limit are all 
planar, and thus do not contain crossed color lines. Moreover, the (single) fermion 
line must always be at the boundary spanned by the planar graph. Consider now a 
typical planar Feynman diagram in the perturbative expansion of ( j2.2| ). Tracking the 
color indices through the diagram, we see that the rectangular off-diagonal blocks of 



272|) vanish identically 

G?(z)=& a (z) = 0, (2.6) 
while the diagonal blocks are proportional to unit matrices, 

G){z) = g M (z) 5), G%z) = g N (z) 5$ . (2.7) 

Upon comparison with fll.8| ) we clearly have 

g N (z) = zG N (z 2 ) and 9m(z) = zG M {z 2 ) ■ (2.8) 



Figure (2) shows some of the diagrams which contribute to G l j(z) to leading order in 
large M, N. 



o = 





a a k k a a 



+ 



Fig. (2) - The first few planar diagrams that contribute to Gj to leading order in the large N, M 
approximation. 



A simple direct calculation shows that the off diagonal blocks of l/(z — H) are 
odd in C, and thus do not contribute to Q2.2p , independently of the perturbation 
expansion. This conclusion clearly remains valid if we generalized ( |2.3| ) into any 
other probability distribution which is even in H, for example, a distribution of the 



form P(C) = ^exp 
defined as usual by 



TrV (C^C) . The self-energy corrections £jv(<z) and Sm(z) are 



9n{z) 



1 



9m{z) 



(2.9) 



: - Y> N {z) ' z- Em{z) 

and correspond to the sum over all one quark irreducible graphs contributing to ( |2.2| ) . 
namely, the amputated quark propagator. For the Gaussian distribution ( [2.3|) , the 
propagators gN,9M and self-energies £tv,£m are related by the simple Schwinger- 
Dyson identities which we display diagrammatically on Fig. (3). 
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Fig. (3) - The two Schwinger-Dyson identities. 



We thus have from Figs. (1) and (3) 



EivO) 



m 2 \fNM 



^ m 2 



9m (z) 



and similarly, 



m 2 \/r 



9n{z) 



(2.10) 



(2.11) 



We substitute the last two equations into (|2.9| ) and obtain the two coupled equations 



9n{z) 



z - ^9m(z) 



-i 



(2.12) 



, 9m(z) =[z- ^t^9n(z)\ 

for cjjy and gu- These two equations clearly transform one into the other under r — > -, 
which interchanges and gu (this is as it should be if the normalization factor in 
(|2.3| ) is symmetrical under M <-> iV, that is, if m 2 is independent of r.) By definition, 
both propagators behave as - in the asymptotic region z — > oo. This asymptotic 
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behavior picks up the physical solution of the quadratic equations for and qm, and 
we find 



9n(z) 



9m{z) 



(a - by 



z 2 -ab- yj(z 2 - a 2 ) {z 2 - b 2 ) 



where 



(a + b) 2 z I 



1 / i _i- 
— (r 4 + r * 
m 



z 2 + ab- yj(z 2 - a 2 ) (z 2 - b 2 ) 



7 1,1 

o = — (r« 
m 



f 4 



(2.13) 



(2.14) 



Note that b measures of the deviation of rectangles from squares: for r = 1, b vanishes. 
The Green's function ( |2.1| ) is thus given by 

1 9N(z)+rg M (z) 



G{z) 



N + M 



N+M 



1 + r 



(2.15) 



where we used (|2.2|) and ( p.7| ). Finally, upon substituting ( |2.13j ) into ( |2 . 1 5|j we find 
that the averaged Green's function of H is 

2 1 



G(z) 



+ b 2 



b 2 ){z 2 



(2.16) 



Let us inspect now some of the features of (|2.16j ). As we discussed in the intro- 
duction, H has M — N "kinematical" zero eigenvalues, regardless of any ensemble 
averaging. In contrast, C^C on the average does not have any zero eigenvalues as 
we have already discussed. Thus, by definition, Gn,m(z) has a simple pole at z = 



with residue ff^ = which is the first term on the right side of (|1.9|). As we 



can see from ( p,14|) , our expression (|2.16|) clearly satisfies this condition provided 
yj [z 2 — b 2 )(z 2 — a 2 ) — ► —ab as z — > 0. This sign of the square root corresponds to 
drawing all four cuts associated with the square root to the left of the branch point 
out of which they emanate. In addition, (|2.14|) is consistent with the required asymp- 
totic behavior - of ( |2.16| ) as z — > oo. The averaged eigenvalue density of H is the 
discontinuity in (|2.16|) as we cross the real axis, except for the origin, which contains 
the "kinematical" zero eigenvalues of H. It is therefore given by 



P(A) 



r + 1 



6(X) 



2 9 [(a 2 - A 2 ) (A 2 - b 2 )} 



7T 



A| 



a 2 + b 2 



[a 2 - A 2 ) (A 2 - b 2 ) 



(2.17) 
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The Green's function (|2.16| ) corresponds to the Hamiltonian ( |1 . 1| ) . With only 
little more effort it is possible to generalize our discussion to calculate the Green's 
function of the Hamiltonian 

H ={c -!)• < 2 ' i8 > 

where e is a fixed "energy" . The off-diagonal blocks fluctuate as before. This Hamilto- 
nian may describe, for example, tunneling between two energy levels with degeneracies 
TV and M that are seperated by an energy difference 2e. For such a Hamiltonian (|1.5|) 
is modified into 

G NM (z) = ^G N (w) + r ^ Z ~f ) G M {w) (2.19) 



r + 1 



r + 1 



where 



2 2 

w = z — e 



such that the identifications (|2.8| ) become 



g N (z) = (z + e)G N (w) 



and 



9m{z) = {z- e)G M {w) . 



(2.20) 



The bare quark propagator in Fig. (1) is split into two pieces, namely, for quarks 
carrying a U (N) color index and for quarks carrying a U(M) index. The defini- 
tions in ( p.9| ) change accordingly into 

1 , , 1 



9n{z) 



z — e — T,m(z 



9m(z) 



z + e- T, M (z) 



(2.21) 



The Schwinger- Dyson identities ( |2.10| ) and ( ^.11| ) are unchanged. We note that the 



set of equations (|2.10|) , (|2.11|) and (|2.21|) are now invariant under 



and 



-e, 



which permutes the two energy levels in H and thus interchanges and gu- With 
these observation it is straightforward to see that ( |2.13[ ) becomes 

2 1 



9n{z) 



(a-b) 2 (z-e) 



w 



ab — \J {w — a?) (w — b 2 ) 



9m{z) 



(a + b) 2 (z + e) 



w + ab — \J(w — a 2 ) (w — b 2 ) 



(2.22) 
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with the same a and b as before. Thus, Q2.15 ) and ( |2.16| ) finally become 



G{z) 



w 



(w 



w 



b 2 ) 



abe 



(2.23) 



(a 2 + b 2 ) w 

which is manifestly invariant under a permutation of the two energy levels of H. 
Note that the matrix ( [2.1§| ) has precisely M — N "kinematical" (i.e., independently 
of the ensemble for C) eigenvectors which correspond to eigenvalue — e. This means 



that ( |2.23j ) has a simple pole at z 



-e with residue and no pole at z 



+e. This property holds provided \/(ir — u- 
already encountered in our analysis of the e 
corresponding to ( [2.23|) is therefore 

• - 1 



(w — b 2 ) — > — ab as w — > which we 
= case. The eigenvalue distribution 



P(A) 



r + 1 



5(\ + e) + 



9 [(a 2 + e 2 - A 2 ) (A 2 - b 2 - e 2 )] 2 |A| 

7i (a 2 + b 2 ) (A 2 - e 2 ) 

In the limit (me) — > oo, the randomness in 



a 2 + e 2 -A 2 ) (A 2 -6 2 -e 2 )(2.24) 



is suppressed, and 



should 



reproduce the eigenvalue distribution of the deterministic part of ( |2.18|) . This is 
indeed the case. In this limit we have -, - — > so both lobes in ( p.24| ) shrink. Each 
lobe contains N eigenvalues, whose number is preserved as the lobes shrink to zero 
width, and thus produce 5 function spikes of strength M ^_ N each. The right lobe 
produces in this way a spike at A = e, while the left lobe coalesces with the already 
existing 5(X + e) spike in (|2.24j) which contains M — N eigenvalues, and thus produces 
a spike containing M eigenvalues. 
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3 Blocks with independent matrix elements and 
their renormalization group analysis 



It is rather difficult to apply the diagrammatic method and sum all the planar dia- 
grams that contribute to G(z) in case of non-Gaussian ensembles. For such ensembles 
that are invariant under the action of some symmetry group one may invoke other 
methods. However, these methods are inapplicable to ensembles lacking the action 
of a symmetry group. 

A class of block structure random matrix models that is not unitary (or orthogo- 
nal) invariant involves matrix blocks C in which each matrix element Cj a is randomly 
distributed independently of the others, with the same distribution for all matrix 
elements. We normalize the matrix elements Ci a symmetrically with respect to M 
and N, such that the two-point correlator 

a 2 

(CiaCjp) = ^——Sijdap (3.1) 

of this probability distribution would coincide with the two point correlator of the 
Gaussian distribution dfj,(C) ~ exp — ^^ TiC^C . For notational simplicity we 
replaced here the m~ 2 in fl2.3|) by a 2 . For this class of matrix models the method of 
orthogonal polynomials is not directly applicable, and alternative methods should be 
sought for. 

For concreteness as well as for simplicity, we consider below the probability dis- 
tribution in which Cj Q may take one of the two values 

± °—ni ( 3 - 2 ) 

with equal probability, where a is a finite number. However, it will be clear from the 
discussion below, that our conclusions are not limited to this particular distribution. 
In order to keep our formulas generic, we therefore treat the Cj Q as complex numbers, 
as long as we do not utilize (|3.2|) explicitly. 



For this ensemble |C ia ,| 2 = deterministically, and thus in particular the 
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diagonal matrix elements of C^C and CC^ do not fluctuate and are given by 
(C ] C) aa = a 2 ^ and (CC% = a 2 /^. 



(3.3) 



We use this convenient property of (|3.2|) in our calculations below. This, however, 



is done with no loss of generality, because in the generic case the non-fluctuating 



quantities in (|3.3| ) should simply be replaced by their averages, which are given on 
the right hand side of (p.3|). 

The random matrix distribution we consider is a generalization of the very first 
model studied by Wigner pUl into random matrices with block structure. Indeed, 
Wigner originally considered large N x N random Hermitean matrices 0, whose ele- 
ments 4>ij = 0*j were either or with equal probability. This matrix model 
follows a semi-circle law for the density of eigenvalues. This semi-circular profile of 
the eigenvalue density was rederived recently JL5J using a large N "renormalization 



group" inspired approach In what follows, we apply the same method to find 
the eigenvalue density p(A) of the random block matrices with independent entries 
introduced above. 

We are interested in 

G(z)= lim G NM {z) (3.4) 

from which p(A) = ^ImG(A — ie) may be extracted immediately. We start our "renor- 
malization group" calculation by trying to relate Gn+i,m(z) to Gn,m(z). To this end 
we consider the M x (N + 1) block 

C = (C, v) 

where v is an M dimensional vector. By definition, each element of C may take now 
one of the two values =fc 1/4 with equal probability. A comparison with the 

original N x M block suggests then that we may draw the C' ia from ( |3.2| ) provided 
we rescale the a parameter in that equation into 
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The non-fluctuating norm squared of v is then given by 



(C' ] C'\ 



N+l,N+l 



(3.6) 



Following we obtain^ after some straightforward algebra 

1 



(N + l)G N+1 (w 
= Tr 



Tr 



(iV+l) 



W 



(JV) 



w - etc - ct^MLc w -v*v -v^C- 

W—V'V ' 



-etc 



(3.7) 



where u> 



We now average over the distribution governing C", keeping terms up 



to O(N ). We first average over the components of v. 

Expanding the two fractions on the right hand side of ( |3.7|) into geometric series 
we see that we have to average products of the form v*Vjvlvi ■ ■ ■ v*v q , which contract 
against products of elements of matrices independent of the V{ . Clearly, 

J 2 



< v*vj >=< C^ +1 C' j: 



N+l 



> 



a 



simply produces a single trace, multiplied by J NM ■ The next non-vanishing correlator 



is 



a 



i A 



(31 



where dijki = 1 when % = j = k = I and otherwise. The last piece in (|3.8|) 
is by definition the fourth order cumulant of the distribution, added to the usual 
pairs of Wick contractions. The correlator ( |3.8| ) then contracts against two matrices 
in the geometrical series mentioned above, producing a term of order f -^=L^ . In 
the large N, M limit, the dominating term in this contraction is the term with the 
maximal independent index summations, which amounts here to two traces. These 
two traces are produced here only by a single pair of Wick contractions. The other 
pair of Wick contractions (which produces only a single trace) as well as the fourth 
order cumulant are therefore negligible in the large N, M limit, and may be discarded 
to leading order. This structure persists for correlators of higher order. The 2n 

2 The analogue of v^v in |L5j was a quantity of 0{jf) which was therefore neglected in the N — > oo 
limit. Here v^v is a finite number and must be therefore retained. 
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point correlator produces in the geometric series a term proportional to J NM ) • m 
that term, a unique string of n Wick contractions produces the maximal number n 
of traces, and therefore dominates the large N, M limit. At this point it becomes 
clear why our calculation, and therefore, the results it leads to are insensitive to the 
details of the distribution of the Ci a . Clearly, only Wick contractions dominate these 
averages in this limit, and thus only the two point function (|3.1|) of the distribution 
matters. This insensitivity to the detials of the distribution was checked explicitly 
in [15 , where various distributions led to the same result. The leading terms in the 



geometric series may be resummed, and one finds the v average 

1 1 



(Tr 



(iV+l) 



W 



etc 



7>« = Tr 



(N) 



w 



etc 



9 , 

ow 



n' 2 a' 2 
w + (N — M)^= - -== w Tr. m 



w 



-etc 



(3.9) 



Invoking large iV factorisation, we can average over the remaining block C immedi- 
ately, by replacing inside the logarithm in (|3.9|) by its average. Thus, 



(N + l)G N+1 (w) = NG N (w) + ^log 



w + 



1-rW 



/ 2 



a 



I 2 



w G N (w) 



. (3.10) 



Remarkably, in the large N, M limit, the v average of Gn + i(w) involves only Gn(w), 
and thus ( |3.10|) is indeed a local (along the N axis) recursion relation involving only 
Gn type Green's functions. This means that the large N "renormalization group" 
recursions for the full Green's function Gn,m close among themselves as we now show. 

Combining ( |1.9|) and ( 3.10 ), we obtain the recursion relation for the complete 
averaged Green's function (|1.4|) 



(N + M+ l)G N+1 , M (z, a') - (N + M)G N , M (z, a') 



oz 



z + 



r)a 



I 2 



[r + l)a 



! 2 



-G N)M {z, a') 



(3.11) 



2ztJ? 2y/r 

where we have displayed the explicit a' parameter associated with the larger C block. 
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In the large N limit (|3.5|) becomes a' — a— ^?+- • •• In this limit, the only possible 



explicit N, M dependence in Gn,m is through the finite ratio r = Therefore we 



may write the left hand side of ( |3.11| ) as 



[(N + M + l)G N+1>M (z, a') -(N + M)G NM (z, a)] + (N + M) [G N , M (z, a) - G NM (z, a')] 



d (N + M) d 

— [(N + M)G NjM (z, a)] + 4^- L ^G NjM (z, a) 



d r + 1 d 
= G N>M (z, a) - r(r + l)—G NM (z, a) H —<r—G NtM (z, a) . 

To leading order in we may drop the N, M indices of the Green's function, replacing 
it by its asymptotic limit (|3.4j ), and replace a' by o inside the logarithm in (|3.11|) . 



The recursion relation (3.11) thus becomes a partial differential equation 



G(z, a) - r(r + l)^G{z, a) + T -^a^G{z, a) 



(3.12) 



d_ 

dz 



log 



2zy/r 



2^ 



(3.13) 



It is easy to see from (1.4) and ( |3.4|) that G(z, a) satisfies the simple scaling rule 

1 „ fz 



which implies that 



G(z,<r) = -G(-,l) 



d d 
a—G(z, o) = -z—G(z, a) - G(z, a) . 
da oz 



(3.14) 



(3.15) 



Thus, using (PUD to eliminate cr£G(z, a) from (|3~lp we arrive at the final form of 
our differential equation for G(z, a), namely, 

^-G(z, a) - r(r + l)^G(z, a) - ^^(z, a) = 



7T lo g 

oz 



z + 



(l-r)a 2 (r + l)cr s 
2zy/r 2Vr 



G{z,a) 



(3.16) 
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This equation tells us how a change in z can be compensated by a change in the 
rect angularity r. 

As a consistency check of our results we can repeat the recursive procedure dis- 
cussed above, but instead of adding an M dimensional column vector to C, we add 
to it an N dimensional row vector u, creating an (M + 1) x N block C" 

C 



C" 



u 



The recursion relation in this case connects, in the large N,M limit, Gm+i{w) and 

Gm{w), and therefore relates Gn,m+i{ z ) to Gn,m(z)- Thus, we simply interchange 

N <-> M in all steps of our calculation above, namely, r <-> -. The differential equation 

for G(z, a) we derived from this recursion reads 

\ c) T ~\~ 1 d 

~^-G(z, a) + (r+ 1)-G(s, a) - _*_ G (z, a) = 



7T lo S 

(72; 



r c 



(r + 1W 



2 — 



G(z,a) 



(3.17) 



2^ 2v/r 

which is indeed the transform of ( p,16|) under r — > i. 

The fact that G(z, cr) satisfies both ( |3.16| ) and its tranform under r — > - means 
that G(z,a,r) = G(z,a, -). This r inversion symmetry of G should be anticipated 
from our N, M symmetric definition of the probability distribution ( |3.1| ) in the first 
place. An important consequence of this r inversion symmetry is that J^G vanishes^at 
r = 1. Thus, at the point r = 1, i.e., for Hamiltonians made of square blocks, (|3.16|) 
reduces to the differential equation 



d 

G(z,a) - z—G(z,a) 



2^1og [z-a 2 G(z,a) 



(3.18) 



previously derived in [TjJ, as expected. 

As was stated at the beginning of this section, only the two point correlator (|3.1|) 
of the random matrix distribution was relevant in the derivation of ( 3.16|) . Hence, the 

r ~ 2 -Q^rrf ( r_1 )i and therefore /'(l) = 



3 This is simply because if f(r) — f (r 1 ) , then £-f(r) 
/'(l) = o. 
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Green's function G(z) of any distribution obeying fl3.1| ) is a solution of ( |3.16| ). We 
have thus shown that for the Wigner ensemble G(z) and the density of eigenvalues 
are universal. In particular, the complex Hermitean distribution ( |2.3j ) as well as 
the real symmetric distribution ( |4.16| ) of the previous section respect (|3.1|) upon the 
identification m 2 = a~ 2 . Thus, their Green's function ( |2.16| ) must be a solution of 
(|3.16|) . A simple check verifies that this is indeed the case. Therefore, the density of 
eigenvalues p(A) = ^ImG(A — ie) is given by ( [2.171) . 

As yet another example of the usefulness of the large N renormalization group we 
use it to prove the central limit theorem in Appendix A. 



By a simple power counting argument (see Section 2 of [13], and also |T7|]) it is 
straightforward to extend the diagrammatic method of the previous section to treat 
the probability distribution considered in this section as well. 
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4 Dyson gas approach 



In this section we present the Dyson gas approach to study the eigenvalue distribution 
of matrices made of rectangular blocks. After completing our calculations we realized 
that our results were previously obtained by Periwal et al. in |Tj|. We assume that 
the M x N rectangular blocks Ci a of the Hamiltonian H in (|1 . 1|) admit the action 
of some symmetry group. Here we focus on blocks with complex entries, but we will 
state some results concerning blocks with real entries in the end. The complex blocks 
are endowed with the natural U(M) x U (N) action 

C -> VCU , V G U(M) , U G U(N) . (4.1) 

One can use this action to bring C to the form 

C=( ^ J (4.2) 

where Ajy is a real diagonal N x N matrix diag(Ai, • - • , Ajv). Therefore, the Hermi- 



tian matrix H in (|1 . 1] ) is a generator of the symmetric space U(M + N)/U(M) ® 
U(N). From these considerations it is clear that C^C may be diagonalized into 
diag(Af , • • ■ , A^v) and CC^ into the same form, but with additional M — N zeros, in 
accordance with (|1.6Q . The probability distribution has to be invariant under (fOj). 
Here we consider distributions of the form 

1 



P(C) = - exp -^MN Tr V{C ] C) (4.3) 

where V is a polynomial and Z is the partition function of these matrices. 

We are interested only in averages of quantities that are invariant under ( f4.1|) . 
We thus transform from the Cartesian coordinates C ia to polar coordinates V^-, U a p 
and X a . Integrations over the unitary groups are irrelevant in calculating averages of 
invariant quantities, which involve only the eigenvalues s a = A^ of C^C. 



The partition function for these eigenvalues then reads ]T8 

N °? N 



Z=l[ [ds a exp[-VNMV(s a )}l[sf- N J] & - s s ) 2 . (4.4) 

a=l q (3=1 1<7<<5<7V 
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The last two products constitute the Jacobian associated with polar coordinates. In 
particular, n( s 7 — s <5) 2 is the familiar Vandermonde determinant. The other product 
is a feature peculiar to non-square blocks. As a trivial check of the validity of (O), 
note that the integration measures in (|2.4| ) and (|4.4j) have the same scaling dimension 
under C -> £C, £ > 0. 

Following Dyson, we observe that ( |4.4| ) may be interpreted as the partition func- 
tion for a one dimensional gas of particles whose coordinates are given by the eigen- 



values s a . The integrand in (|4.4|) may be expressed as exp [— yNM S] where 

1 



£ = E 



N r _ 1 

V(s a ) ■=- log s c 



a=l 



y rr. E log (*« " *tf ( 4 -5) 

iV V l<a<l3<N 

is the energy functional of the Dyson gas. In the large N, M limit (|4.4|) is governed by 
the saddle point of (|4.5|) , namely, by a C^C eigenvalue distribution {s a } that satisfies 

dS _ , r-1 1 2 1 

ds a a yjr s a N^/r ^ s a - s p ' 

Here the prime over the sum symbol indicates that j3 = a is excluded from the sum. 

We now turn our attention to the average eigenvalue density of H, which we may 
readily deduce^]]] from the averaged Green's function Gn,m{ z ) m $01) • The s a are 
eigenvalues of C^C. We thus calculate first Gjq{z 2 ), which according to flL~8| ), is given 
by 

1 N 1 

Gn(w) = -Y,( )• (4-7) 



Here the angular brackets denote averaging with respect to (|4.4j) . By definition, 
Gn(w) behaves asymptotically as 

G N (w) — ► -. (4.8) 

w— >oo yj 

It is clear from (|4.7|) that for s > 0, e — ► 0+ we have 



1 * 1 CT w 

Gn(s - ie) = -P.P. E<— -> + T7 - (4-9) 

iV a=l S S " iV af=l 
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where P.P. stands for the principal value. Therefore, the average eigenvalue density 
of C^C is given by - Im Gn{s — it). In the large N, M limit, the real part of ( {Op is 
fixed by (|4.6|), namely, 



Re G 



N(S - 16 



V'(s) 



(4.10) 



The potential V(s) in (|4.3|) clearly has at least one minimum for s > 0, and will 
therefore cause the eigenvalues to coalesce into a single finite band or more along the 
real positive axis. Moreover, the log s term in (|4.5|) clearly implies that the {s a } are 
repelled from the origin. We thus anticipate that the lowest band will be located at 
a finite distance from the origin s = 0. 

At this point we depart from discussing the general distribution and assume for 
simplicity that the probability distribution is given by the Gaussian distribution (|2.3|) 
with 

V(s)=m 2 s. (4.11) 

In this case we expect the {s a } to be contained in the single finite segment < b 2 < 
s < a 2 , with a > b > yet to be determined.^ This means that Gn(w) should have 
a cut connecting b 2 and a 2 . This conclusion, together with Q4.1U| ) imply that Gn{w) 
must be of the form 



Gn(w) 



r m 



+ F(w) ^J(w -b 2 )(w - a 2 ) 



where F(w) is analytic in the w plane (with the origin excluded.) The asymptotic 
behavior (M.8I) then fixes 



and thus, 

G N (w) 



F(w) 



1 

2^ 



r m 



2w 



a 2 + b 2 



1 



77T 



(4.12) 



r m 2 w — r + 1 — yfr m 2 J \w — b 2 )(w — a 2 ) 



(4.13) 



4 We find below, of course, that a and b coincide with the expressions in ( |2.14[ ). 
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The eigenvalue distribution of C^C is therefore 

p(s) = -lmG N (s- te) = J(s-b 2 )(a 2 -s) (4.14) 

7T Z7TS v 

for b 2 < s < a 2 , and zero elsewhere. 

We now substitute G^(z 2 ) from ( |4.13| ) into (|1.9|) to obtain an expression for 
Gn,m( z )- As we discussed in the introduction and in section 2, Gn,m( z ) has a simple 
pole at z = with residue ffqr^ = ^Ej, which is the first term on the right side of 
(|1.9Q . We thus conclude from ( |1.9|) that wGn(w) must vanish at w — 0, which in 
turn implies a second condition^ on a, b, namely, 

ab=-^=. (4.15) 

We are now able to fix a and b from ( [4.12|) and (|4.15| ) and find that they are given 
by ( |2.14| ). We thus find that Gn,m(z) coincides with ( |2.16| ) and that the averaged 
eigenvalue density of H is the expression in (|2.17|) . 

We close this section by sketching the similar analysis of Gaussian random Hamil- 
tonians made of real M x N blocks C. We parametrize the Gaussian real orthogonal 
ensemble by 



1 r m 



2 



P(C) = -exp [-— VNM Tr C T C] (4.16) 

Zj Z 



with the partition function 

M N 



Z= f ]J IT d C ia exp [-— y/NM Tr C T C] . (4.17) 
The two point correlator associated with ( |4.16|) is clearly 

(CiaCtf) = Wop- (4-18) 

Note that (|2.3| ) and ( 4.16| ) are conventionally parametrized in such a way that (2.5) 
and ( |4.18| ) coincide. 



5 Note that the Riemann sheet of the square root in (4.13) is such that y/{0 — b 2 )(0 — a 2 ) = —ab, 
as we already observed in section 2. 
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The partition function for the corresponding Dyson gas reads [18 

N 00 

n 

a=l q ~ /3=1 !<7<<5<A r 



z= n ds « ex p t-^viVM m 2 sj n s p 2 n 1% - . (4.19) 

a=l q ^ /3=1 1<7<5<7V 



As before, the last two products constitute the Jacobian associated with polar coor- 
dinates. The energy functional £ of the Dyson gas is now 

Thus, in the large A/", M limit, ( [4.20 ) becomes precisely one half of the corresponding 
expression (|4.5f) for complex Hermitian matrices, and our discussion following (|4.6|) 
through (|2.17| ) remains intact. 



23 



5 Kazakov's method extended to rectangular com- 
plex matrices 



5.1 Contour integral 

Gaussian matrix ensembles may be studied in many ways. Several years ago, 



Kazakov introduced a method for treating the usual Gaussian ensemble of ran- 
dom Hermitian matrices, which was later extended and applied to a study of random 
Hermitian matrices made of square blocks0. Here we generalize it to random Her- 
mitian matrices made of rectangular blocks. It consists of adding to the probability 
distribution a matrix source, which will be set to zero at the end of the calculation, 
leaving us with a simple integral representation for finite N. As we will see, one 
cannot let the source go to zero before one reaches the final step. We modify the 
probability distribution (|2.3|) of the matrixF] C^C by adding a source A, an N x N 
Hermitian matrix with eigenvalues (ai, • • • , a at) : 



P A (C) = — exp(-v / MiVTrC t C - VmNTtAC^C). (5.1) 
Za 

Next we introduce the Fourier transform of the average resolvent with this modified 
distribution: 

U A (t) = <^Tre^ tc ), (5.2) 
from which we recover the eigenvalue density 

_ w ^e- Us U (t) = (jj™ (a - &C)) (5.3) 

of C^C, after setting the source A to zero. Without loss of generality we can assume 
that A is a diagonal matrix. Let us now calculate £/^(t). We first integrate over 
the N x N unitary matrix U which diagonalizes C^C . This is done through the 
well-known Itzykson-Zuber integral over the unitary group [22 

r det\e aab ' 3 ] 
J^MVAUBUl) = (5.4) 

6 For notational simplicity we set m 2 — 1 in ( |2.3|) throughout this section 
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where A(A) is the Vandermonde determinant constructed with the eigenvalues of A: 



A{A) = J] (a a - a p ) 

a</3 



(5.5) 



(pi, - • ■ , 6jv) are the eigenvalues of B, and A(B) is the Vandermonde determinant built 
out of them. We are then led to 

| .Y 

Z A A( 



u A (t) 



i i N r 

]y E J dsx---ds N e ltSa A(si, • • • , s N ) 



x 



2V 

\p=i 



M-N 



cxp 



iV 



-VMiV^s 7 (l + 



7=1 



(5.6) 



We now integrate over the s a 's. It is easy to prove (for example, by using the Faddeev- 
Popov method) that 

M-N N 

exp(- ^ s a 6 Q ) 



y dsi ■ • ■ ds N A(si, ■ ■ ■ , s N ) I f[ S/3 
A(6i, • • • , M 



v/3=l 



a=l 



(5.7) 



yJV (nf&«)« 

where Cat is a constant independent of the b a . Note that ( [5.7D is valid also for M = N. 
With the normalization Ua(0) = 1, we could always divide, at any intermediate step 
of the calculation, the expression we obtain for U A {t) by its value at t — 0, and thus 
the overall multiplicative factors in (5.6) and (5.7) are not needed. 
We now apply this identity to the iV terms of ( |5.6j ), with 

it 



bf(t) = VMN(l + a p 



(5i 



and obtain 

u A (t) 



-t N N i i n 

iv a=ip=l L 



a 



n 



CLP Cb^y 



>MN 



{5a,/3 — $0,7) 



it x 

MN Ua >P /3<7 



dp — CL~ 



1 ^ ^ r 1 ~\~ d(X 

ivi l l + a«-- 



n 



MiV 7t^ q 



- Qq Vmn 



(5.9) 
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As a consistency check, note that for M — N, ( |5.9| ) coincides with Eq. (3.9) of J7J. 
This sum over N terms may be conveniently replaced by a contour integral in the 
complex plane: 



U A (t) 



* r du 



t 



1 + u 



M 



N 



2ni \ 1 + u 



n 



u — 



fMN 



U — QU 



(5.10) 



fMN J 7=1 T 

in which the contour encloses all the a 7 's and no other singularity. It is now, and 
only now, possible to let all the a 7 's go to zero. We thus obtain a simple expression 
forf/oW, _ 

i \f9 r du (i - - 



W) 



/A/TV 



A' 



2m 



M 



(5.11) 



(i+«)VMAf y 

Note that this representation of Uo(t) as a contour integral over one single complex 
variable is exact for any finite M, N, including M — N — 1. 



5.2 The density of states 



In the large M,N limit (with finite r = 4^-), for finite t, the integrand in Q5.ll 



becomes e 



and therefore Uo(t) approaches 



Jr f du it 



where we changed u into —it. 



Setting z 



t^- we change variables to 

l—u ° 



U 



2z 



Then the integral of ( |5.12 ) becomes, after an integration by parts, 



TT (+\ J T I dZ dU Jtz 



dz 



-Vr j> ^ u{z) e 



itz 



2m 



= ^ — J {a? -x)(x- 6 2 ) e 4te 
2vr 7&2 x v v ; v y 



1 



>/r + 



4z 



(5.12) 



(5.13) 
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(5.14) 
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where a and b are given in ( 2.14|) . Therefore, we have from ( |5.3| 

dt 



2,71 



P(s) = 

rz .Fui - 1\ r« -hi\ 

(5.15) 

for b 2 < s < a 2 , and zero elsewhere. This expression coincides with ( f4.14|) as expected. 
We observe from ( |1.9|) that p(A) and p(s) = p(A 2 ) are related by 



P(A) = ^— T<KA) + ^Lp(A 2 ) . 
r + 1 r + 1 

Substituting ( |5.15| ) into ( [5.1 6| ) we obtain ( |2.17| ) once again, as we should. 



(5.16) 



5.3 The edges of the eigenvalue distribution 

It is easy to apply this same method for studying the cross-over at the edges of 
the eigenvalue distributions ( |2.17| ) or ( |5.15| ), namely, in the vicinity of the end points 
s = a 2 and s = b 2 . To this end, we observe from (|5.3| ) and (|5.11|) that 

N 

(5.17) 



ds 



dt 



-its 



2vr 



du 1 



Nuy/r 



2iti A _ i±/r n M 

M(l+u) 



where the purpose of the s derivative is to get rid of the simple pole at t = in ( |5.11| ). 
By changing t to yMN t and then t to t + iu, as well as u to — iu, we obtain the 
factorized expression 



dp(s) 



-iM 



°° dt 



-WMN ts 



t 



N 



du 



JVMN us 



[u + i 



,Af ■ 



(5.18) 



ds \J-oo2n~ (t + i) M ) \J 2wi u 1 

The advantage of (|5.18|) is that it is relatively easy to study its large N, M behavior 
by saddle point techniques. We observe that the t integral may be written as 

/+oo 
dte -VNM S eff (5 lg) 

-co 

where S e ff is given by 

S eff = i s t + log (t + i) - -^= log t . ( 5.20 ) 
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Similarly, the integrand of the u integration is e^ NM Seff . Thus, the large N,M be- 
havior of (|5.17|) is determined by the saddle points of a single function S e ff. Consider 
the t integral ( |5.19| ) first. It has two saddle points t c at 



s-ab±^/(s-a 2 ) (s-b 2 ) 
2is 



(5.21) 



where a, b are given in Q2.14|) . The interesting situation occurs when these two saddle 
points become degenerate, namely at the endpoints s = A 2 = a 2 and s = X 2 = b 2 . We 
thus investigate ( J5.17|) at the vicinity of these points, by focusing on these regions. Let 
us consider the neighborhood of A = a first (the cross over behavior around A = —a 
is simply the mirror image thereof.) We introduce the scaled variables 



A 
t 



a + N~ a x, 
1 



+ N~ P T 



with a,j3 to be determined, and expand S e tt up to r 3 . This leads to 



S eff {t) = S* + 2 r-4 N~ a x 

+ - a 4 r 3 N~ 3 P + 2ia N'^rx + ■■■ 
3 



(5.22) 



(5.23) 



where S* is the value of S e ff at the critical point, and the ellipsis stand for terms of 
0(N~ 2a ). We thus find that there is a large N, finite x limit, provided we fix the two 
unknown exponents a and (5 to 



2 

3' 



B= l - 

H 3 



(5.24) 



We repeat this for the u-integral of ( ^.18|) . We then find that the leading terms of 
(|5.23|) of order 1, as well as the term 2xN~ 2 ^ 3 , cancel with terms of opposite signs in 
the ^-integral. Thus we obtain the following equation for the density of state near 
the critical value s = a 2 or A = ±a, 



dp(\ 2 
d\ 2 



-Ms 



Ai 



2 [ -VjV§(A=Fa) 
a) 



(5.25) 
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where the Airy function Ai(z) is defined as 

Ai\{3ay 1/3 x] = ± — '- — / cos(at 3 + xt)dt. (5.26) 

7T JO 

The Airy function in ( |5.25| ) is smoothly decreasing for |A| > a but it oscillates for 
|A| < a. 

Investigation of the behavior of ( 5.18Q near the other critical points A = ±6 pro- 



ceeds similarly. Concentrating on A = b we introduce the scaling variables 

A = b + N- a x, 

t = —L—+N-PT (5.27) 
l r — i 



and find that there is a large N, finite x limit, provided we fix the two unknown 
exponents a and (3 to the same values as before. Thus, the crossover behavior of the 
density of states arround A = ±6 is governed by the Airy function as well, for any 
r > 1. 

A new phenomenon appears, however, if we also take the limit r — > 1. It is easy 
to see, by rescaling r in the expansion of S e ff into 

T=(v^-l)r, (5-28) 

that the Airy function behavior of 9 ^ 2 ~ near A = ±b breaks down as r — > 1. Indeed, 
from previous work [l], |2|, [3], |], |5], [| 0, § we know that the oscillations near the origin 
in the density of the eigenvalues of matrices built out of square blocks (r = 1) are 
governed by the Bessel function and not by the Airy function. 



ACKNOWLEDGEMENTS 

This work is supported in part by the National Science Foundation under Grant 
No. PHY89-04035. 



29 



Appendix : The Central Limit Theorem - A Renormalization Group 
Proof 



As a simple, but perhaps amusing exercise we use the large A" renormalization 
group discussed in Section 3 to prove the celebrated central limit theorem of Gauss. 

Consider a set of A" independent random variables {xi, X2, ■ ■ ■ , xn} which are 
distributed according to some distribution function 

N 

Qn(xi, • • • ,x N ) = Y[Q(xi) . (A.l) 

i=i 

In order to be consistent with our normalization conventions in Section 3, we normal- 
ize this distribution function such that 



where j3 > is yet to be determined. Thus, a typical term drawn from Qn{x) is of 
the order a N' 13 . We wish to calculate the distribution function of the sum of these 
random numbers, namely, the quantity 



where (-)n denotes averaging with respect to Qn(x). In principle, Pn depends upon 
all the cumulants of Qn(x), but we expect that the large A^ limit of Pjv will depend 
only upon a. Following our discussion in Section 3, we now consider a set of A^ + 1 
random variables whose distribution function Qn+i(x) is normalized such that 





AT2/3 °H 



(A.2) 




(A.3) 



{xi) = 0, 



(XiXj) 



(A.4) 



Then, 





+ 




(A.5) 
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where we used ( |A.4j ) . The ellipsis stand for cumulants of order higher than two, which 
are clearly suppressed by powers of N~@, and we neglect them henceforth. Comparing 
and (|A.4p we also see that 



5 (s-^x^J ) n+i = P N {s,a') 



with 



a 



N 



a = (l 



p_, 

.N+U K ~ N' 

We now use ( A.6 ) and ( |A.7| ) to rewrite ( |A.5| ) as 



a 



P N+1 (s,a) 



1 



P d 



o -7T- + 



a 2 d 2 



P N (s,cr) 



(A.6) 



(A.7) 



(A.8) 



N da 2N 2 P ds 2 

where we neglected terms of 0( N2 ], +1 ). We observe from (|A.8| ) that variations of a 
are as important as variations of s in the large iV limit only if 

1 



P 



(A.9) 



which fixes p. We thus conclude that 



N 



dP N a 



N 



d 2 







a 



ds 2 da 



P N {s,a) 



The left hand side of ( |A.9|) must vanish if Pjv has a large limit 



P N (s,a) — > P(s,a) , 



and thus 



d 2 d 



a 



P(s,a) = 0. 



(A.10) 



(A.ll) 



(A.12) 



ds 2 da 

A simple scaling argument, similar to the one invoked in Section 3, leads to the 
relation 



P(s,a) = -P(-,1) (A.13) 
a a 



which implies that 



*° p = -p- s §- p. 

da ds 



(A. 14) 
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Substituting ( |A.14j ) in ( |A.12| ) we finaly obtain the differential equation 







(A.15) 



We solve ( |A.15| ) and find that its normalized solution is the Gaussian distribution 



P(s,a) 



'lit a 



cxp 



2a 2 



(A.16) 



which is the statement of the central limit theorem. The proof of the central limit 
theorem presented here is not any simpler than the conventional proof found in text- 
books. 



The generalization of this proof to the case|23|] of adding a large number N of 
K x K matrices {<f>x, ■ ■ ■ , </>at} is straightforward. In this case s and P(s,a) are 
K x K matrices. We take these matrices to be real (the Hermitian case can be 
treated similarly.) Then (|A.15|) becomes 

( o d 2 



(A.17) 



where /z, v are indices of the K x K matrices (repeated indices are summed over.) 
The normalized solution of ( |A.17p is the Gaussian distribution 

Tr s 2 \ 



P{s, a) = (v 7 ^ K a) K exp 



2A^V 2 



(A.18) 
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